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Simulated tempering (ST) has attracted a great deal of attention in the last years, due to its 
capability to allow systems with complex dynamics to escape from regions separated by large entropic 
barriers. However its performance is strongly dependent on basic ingredients, such as the choice of 
the set of temperatures and their associated weights. Since the weight evaluations are not trivial 
tasks, an alternative approximated approach was proposed by Park and Pande (Phys. Rev. E 76 , 

016703 (2007)) to circumvent this difficulty. Here we present a detailed study about this procedure 
by comparing its performance with exact (free-energy) weights and other methods, its dependence on 
the total replica number R and on the temperature set. The ideas above are analyzed in four distinct 
lattice models presenting strong first-order phase transitions, hence constituting ideal examples in 
which the performance of algorithm is fundamental. In all cases, our results reveal that approximated 
weights work properly in the regime of larger R’s. On the other hand, for sufficiently small R its 
performance is reduced and the systems do not cross properly the free-energy barriers. Finally, 
for estimating reliable temperature sets, we consider a simple protocol proposed at Comp. Phys. 

Comm. 128 , 2046 (2014). 

PACS numbers: 05.10.Ln, 05.70.Fh, 05.50.+q 


I. INTRODUCTION 

Although Monte Carlo method has become probably 
the most common tool for studying phase transitions and 
critical phenomena, in practice its usage is not so simple, 
whenever standard algorithms (e.g. Metropolis) are used. 
Despite the simplicity and generality, they lead to diffi¬ 
culties close to the emergence of phase transitions. For 
instance, alternative procedures are typically required, 
specially in the case of systems with microscopic config¬ 
urations separated by valleys and hills in the free-energy 
landscape |l-4j- Cluster algorithms [H, [g], multicanon- 
ical Wang-Landau Q and tempering methods are 
some examples of proposals to overcome these difficulties. 
In particular, parallel tempering (PT) [9] and simulated 
tempering (ST) [To| enhanced sampling methods have 
drawn attention due to their generality and simplicity 
when compared with the previous examples. Their basic 
idea consists of using configurations from high tempera¬ 
tures for systems at low temperatures, allowing in prin¬ 
ciple the dynamics to escape from metastable states and 
providing an appropriate visit of the configuration space. 
In particular, distinct aspects of tempering methods have 
been explored in the last years, aiming at better under¬ 
standing of efficiency and pertinence. For instance, the 
role of temperature sets for the PT case was investigated 
in Refs. 11 14j, whilst the importance of non-adjacent 
exchanges was taken into account in Refs. jT~il - iT8^ . In ad- 
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dition the efficiency and comparison between tempering 
methods were considered in Refs. fl9l - l2lj . 

Focusing our attention in the ST we face one of its 
main difficulties, namely the evaluation of the free-energy 
weights, required for an uniform sampling to all tempera¬ 
tures. Despite the development of alternative techniques, 
their applicability for more complex systems still poses a 
hardship. In some cases jUHII: tire accumulation of his¬ 
tograms (of a given quantity) and previous simulations 
are necessary to calculate (or to estimate) the input pa¬ 
rameters that guarantee a sufficient number of visits to 
all temperatures. In such cases, the weights are set ar¬ 
bitrarily but a knowledge of the partition function Zi at 
each temperature is required and a flattening histogram 
based on a random walk in the parameter (temperature 
or energy) space is used to obtain a satisfactory esti¬ 
mation of Z, L . In Refs. [23], [24|, the partition function 
is exactly valued through numerical simulations, taking 
into account its relationship with the largest eigenvalue 
Aof the transfer matrix T. Although the evaluation 
of A(°> is possible for lattice-gas systems, its extension 
for more complex cases (e. g. off-lattice systems) is 
not straightforward. In contrast to previous “exact” ap¬ 
proaches, Park and Pande f2j| proposed an approximated 
tool of estimating weights, based on the average system 
energy. Since the mean energy is easily obtained for any 
system (including lattice and off-lattice models), it con¬ 
stitutes a considerable simplification over the free-energy 
case. Nevertheless, there are some fundamental points 
that need to be understood in order to make it a promptly 
useful method. The first one is how this procedure com¬ 
pares itself with using free-energy weights? The second 
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one is under what conditions does it provide equivalent 
results to those obtained from free-energy weights? An 
additional point is if it is possible to obtain proper tem¬ 
perature set that yields precise results under lower com¬ 
putational cost. To answer the aforementioned points, we 
have analyzed, under the ST with approximated weights, 
four distinct lattice models, namely, Blume-Capel (BC) 
and Bluine-Emery-GrifBths (BEG) [26j], Bell-Lavis (BL) 
and associating lattice gas (ALG) water models 
The former two are interesting tests, due to 
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the existence of very precise results available from clus¬ 
ter algorithms [3lj, Wang-Landau 32| . PT and ST with 


free-energy weights in la min. Therefore, they 


con¬ 


stitute relevant benchmarks for our purposes. The BL 
and ALG are also important examples, taking into ac¬ 
count their more complex phase diagrams, including liq¬ 
uid phases with distinct structures, regions of unusual 
behaviors (density and diffusion anomaly lines) and also 
dynamic transitions [a HU- I n the case of ALG, an extra 
advantage arises, due to the existence of two phase coex¬ 
isting lines, between gas and liquid phases. Hence, the 
ALG works as a double checking of reliability of our pro¬ 
posals. For instance, we focus on the regime of low tem¬ 
peratures, in which strong first-order phase transitions 
separate coexisting phases. Recently a general a ppr oach 
for discontinuous transitions has been proposed [35], [36[, 
in which thermodynamic quantities are described by a 
general function, allowing to achieve all relevant infor¬ 
mations by studying rather small system sizes for some 
control parameters. Thus, its combination with a proper 
usage of ST can provide us a powerful approach to deal 
with discontinuous transitions with rather low computa¬ 
tional cost. 

Henceforth, the analysis of all cases will show that the 
approximated weights work properly (and hence lead to 
correct results) in the regime of large replica numbers 
R for an appropriate choice of temperature sets [33]. On 
the other hand, for sufficiently small R its performance is 
strongly reduced and the system does not visit properly 
the distinct coexisting regions. Finally we extend for the 
approximated weights, a simple protocol for obtaining 
proper temperature sets initially proposed for the free- 
energy weights [33]. 


II. SIMULATED TEMPERING AND 
APPROXIMATED WEIGHTS 

The basic idea of the ST concerns with the fact that 
the system temperature T can assume different values 
between the extreme values Tf and Tr. where R is the 
replica number. The MC simulation is defined as follows: 
In the first part, starting at a given temperature X) within 
the set Tr = {Ti,...,Tr} (in all cases we started from 
Tr), a given site of the lattice is randomly chosen and its 
variable is changed (among all possibilities) according to 
the Metropolis prescription nhn{l, exp(/3,AH)}, where 
A'H(a) denotes the energy difference between the ’’new” 


and ’’old” configurations and fa = 1 /hsT. After re¬ 
peating above dynamics a proper number of realizations 
(here L 2 random choices are considered) in the second 
part the temperature exchange (T) —>• X)) occurs with 
the following probability 


Pi^yj = min{l, exp[(ft - fa)R{a) + (g 0 - &)]}, (1) 


where g\ is the weight associated with the temperature 
Ti and R (a) is the system Hamiltonian. For a uniform 
sampling, the weights should tie proportional to the free- 
energy /,; given by gt = fa fi [25] . Since the evaluation of 
/ is not an easy task, alternative calculations of weights 
have been proposed 121 | - [2-l| . The simplest proposal [25] 
estimates the g 's according to the following approximated 
formula 


9i-9i*{Pj-Pi)<JJj+Ui)l 2, ( 2 ) 

with Ui = {Ri) (i = 1,2 ,R) denoting the average 
system energy at Ti. Thus, from Eq.([2]) the weights 
are estimated from simple and direct standard numerical 
simulations. Here we give a further step by analyzing 
them by inspecting two crucial points: their dependence 
on the replica number R and on the set of temperatures 
Tr. In order to scrutinize them, we compare numerical 
results at the phase coexistence points for distinct R' s, 
with temperature schedules estimated as proposed in Ref. 
[33| and described as follows: starting from a fixed T\ we 
choose the next R— 1 temperatures X 2 < X 3 < ... < Tr 
in such a way that the resulting exchange frequencies 
fi+i t i between any two successive temperatures Tj and 
T i+ i are all equal to some value specified = /. We 
define 1 as the ratio of the number of exchanges be¬ 
tween Ti and T i+ i to the total Monte Carlo steps Nmc- 
Note that from this recipe the highest temperature Tr 
becomes automatically obtained. The efficiency of such 
achieved set Tr is verified by means of standard tests, 
where in the case of first-order transitions, the tunneling 
between the coexisting phases and convergence to the 
steady state starting from a non-typical initial configura¬ 
tion constitute proper efficiency measures. More specif¬ 
ically, the existence of full trapping in a given phase or 
even temperature changes that do not allow the system 
to visit properly the coexisting phases will imply in ther¬ 
modynamic averages marked by no changes or abrupt 
variations (see e.g Figs. 2 and 3 for / = 0.02 and 0.021, 
respectively). Such points can be understood by recall¬ 
ing the ideas from Ref. [35, 36], when the system close to 
the phase coexistence have typical thermodynamic quan¬ 
tities, like energy and order parameter, well described by 
the following general expression 

m a r 

W ( y ) ~ (h + bn ex P h a ny ])/(l + ^2 Cn ex P[ -a ^D> 

n —2 n —2 

(3) 

where for J\f coexisting phases, y denotes the “distance” 
to the coexistence point £* given by y = £ — £*. The 
coefficients c n ’s and 6 „’s are related to derivatives of the 
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free energies /„ of each phase n with respect to param¬ 
eter £ reading d/„/<9£ (36|. In the case of two phase 
coexistence (A f = 2) Eq. © acquires the following 
way W(y) = (&i + b 2 exp[-a 2 y ])/(1 + c 2 exp[-a 2 y]) and 
hence only four parameters are necessary to determine 
the whole function. In other words, according to Eq. 
©, numerical simulations (of a given system) for known 
L and control parameter sets £ = £o (like chemical po¬ 
tential and temperature) will provide a well defined value 
for thermodynamic quantities Wq = W(L,y 0 ), where 
Vo = Co — £*■ Note that at the phase coexistence point 
y = 0, the quantity W reads 


W 0 


bi+TZ= 2 b, 

1 + El 2 Cr 


(4) 


for all Us. Hence, different curves of W should cross at 
the coexistence point. However, as Wo and Wq are ver¬ 
ified only for dynamics that visit properly the distinct 
phases (e.g. one flip algorithms lead to strong hysteresis 
at low T”s and results do not obey Eq. © ). In the 
case of tempering methods, the achievement of results 
not following Eq. © indicates that Tr is not proper. 
Typically, low TVs (including the extreme Tr) provide 
high temperature exchanges, but the system is not able 
to cross energy barriers. Hence / should be decreased, 
in order to raise the temperatures. In contrast, a very 
high 7 ~r (obtained from a very low /) may be respon¬ 
sible for few frequent exchanges and hence poor aver¬ 
ages are obtained. So that / should be raised, in order 
to increase the temperature changes. An intermediate 
optimal frequency f op t is expected to satisfy the above 
points and hence it may lead to frequent tunneling be¬ 
tween the phases and a faster convergence to the steady 
state. Some remarks are needed: Although our present 
recipe provides a reliable temperature set, it does not 
exclude the existence of other optimized choices present¬ 
ing distinct frequencies between adjacent temperatures. 
Also, other efficient sets (ranged from T\ and Tr) can be 
obtained. For simplicity, we have considered such equal 
frequencies proposal. 

In this work, we will test this procedure for distinct 
R's, in order to verify the replica number influence in 
the results. It is worth mentioning that although some 
numerical work is necessary to find / op t, in practice such 
process demands relative short simulations. Hence, the 
search for the corresponding optimal Tr is not a compu¬ 
tationally time consuming procedure. 


III. MODELS 

A. BEG and BC models 


where Oi = 0 , if the site i has null spin and ±1 if i has 
up and down values, respectively. The parameters J and 
K are interaction energies and D denotes the crystalline 
field. The BC model corresponds to the K = K/J = 
0 case. In order to compare with previous results [HI, 
m H3| we also consider the K = 3 case, in which the 
system displays ferromagnetic and paramagnetic phases. 
All phases can be characterized in terms of two order- 
parameters q = (of) and m = (Oi). At low temperatures, 
all transitions are first-order and yield close to the T = 0 
the value D* = z(K + l)/2, where z is the coordination 
number of the lattice. In both model studies, we consider 
a two-dimensional square lattice in which z reads 4. 


B. Bell-Lavis model 

The Bell-Lavis (BL) model is defined on a triangular 
lattice (2 = 6 ) where each site is described by two kinds 
of variables, namely occupational (a) and orientational 
(r) states. The former takes the values a* = 0 or 1, when¬ 
ever the site i is empty or occupied by a water molecule, 
respectively. The latter describes the possibility of form¬ 
ing hydrogen bonds and reads t \ 3 = 0 or t \ 3 = 1 when 
the arm is inert or bonding, respectively. Two molecules 
interact via van der Waals and hydrogen bond energies 
provided they are adjacent and point out their arms to 
each other (r)' J r ) 1 = 1), respectively. The BL model is 
given by the following Hamiltonian 

H = - a, Oj ( e hb rf if + e vdw ) - p ^ cr*, (6) 
<i,j> i 

where e vdw and ehb are the van der Waals and hy¬ 
drogen bonds interaction energies, respectively, and p 
is the chemical potential. We also compared results 
with those obtained from PT and ST with free-energy 
weights [24], [33]. For instance, we consider £ = 0.1 
(where £ = e vdjW /ehb), in which the system presents a 
gas and two distinct liquid phases, named low-density- 
liquid (LDL) and high-density-liquid (HDL) [27],[l8|. In 
the gas and HDL phases, the lattice is empty and it is al¬ 
most filled by molecules respectively. On the other hand, 
the LDL phase presents an intermediate density close to 
p = 2/3 and exhibits an ordering structure (honeycomb 
like geometry), signed by a maximum density of hydrogen 
bonds per molecule. At T = 0 (where T = T/ehb), both 
transitions are first-order and occurs at p* = —3 (l + £)/2 
and p* = — 6 £, respectively. For T^O and distinct £’s 
the former phase transition remains first-order, but the 
latter becomes critical [H[. For £ = 0.1, the second-order 
and first-order lines meet in a tricritical point. 


The BEG model Hamiltonian [2(| is given by 


C. Associating lattice-gas (ALG) model 


n =- yy ( ja i t+ D y2<?h 

<i,j> i 


( 5 ) 


Similarly to the BL, the ALG model is also described 
by an occupation (a) and orientation (r) states. But 
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an important difference from the BL is that an energetic 
punishment exists when a hydrogen bond is not formed. 
More specifically, two first neighbor molecules have an 
interaction energy of — v (—v + 2u) if there is (there is 
not) a hydrogen bond between them. The Hamiltonian 
system reads 

n = 2uJ2 ~vYl ai - 

<i,j> i 

The ALG also presents a gas, LDL and HDL phases, with 
densities p = 3/4 and 1, respectively Another rel¬ 
evant distinction from the BL model is that here both 
gas-LDL and LDL-HDL transitions remain first-order 
for T = T/v ^ 0 and end at respective critical points 
joined by a critical line separating a structured liquid 
from a disordered fluid phase. At T = 0, the discontin¬ 
uous transitions occur at p = p/v = —2 (gas-LDL) and 
p = —6 + 8 u/v (LDL-HDL). As in the BL model, p is a 
proper parameter for the gas-LDL transition, but not for 
other phase transitions. Alternative </>’s are (j> = (4 p — 3) 
or also (f> = (2phb~3), being phb the hydrogen bonds den¬ 
sity, taking the values phb = 3/2 and 1 for LDL and HDL 
phases, respectively. Here, we focus on the temperature 
T = 0.30, in which results from the PT show that the for¬ 
mer and latter transitions take place at p* = —1.9986(2) 
and p* = 2.0000(2), respectively [35:]. 


IV. NUMERICAL RESULTS 

Numerical simulations were carried out for several sys¬ 
tem sizes L and control parameters. In all cases, nu¬ 
merical simulations start from a fully ordered initial con¬ 
figuration and 3.10 6 MC steps have been used to “equi¬ 
librate” the system. After the transient, we have used 
3.10 6 MC steps for evaluating the weights and steady 
analysis start from the largest temperature Tr. Except 
Fig. Ha), all quantities have been evaluated after the 
transient regime. From now on, we will replace all re¬ 
duced quantities (T,H,D,p) for (T, H, D, p). In order 
to achieve a global idea about approximated weights and 
how they compare with free-energy ones, in the first anal¬ 
ysis we show the tunneling between coexisting phases for 
the BEG and BL models, respectively by taking the tem¬ 
perature schedules Tr (for distinct i?’s) obtained from 
free-energy weights [33j. For low temperature values, 
Ti = 0.5 (BEG) and 0.1 (BL), the coexistence points 
yield at D* = 8.0000(1) and p* — —1.6500(1), respec¬ 
tively. In these points, according to Eq. m a correct 
system sampling should give steady equilibrium values 
consistent with go « 2/3 and po ss 1/2 for the BEG 
and BL, respectively. The former can be understood by 
recalling that two ferromagnetic phases (q ~ 1 ) coex¬ 
ist with one paramagnetic phase {q ~ 0). Since at the 
phase coexistence their weights are equal (1/3), we have 
go ~ 2/3 for all L’s. A similar reasoning holds for the lat¬ 
ter model, in which the LDL phase has density p « 2/3 
(with degeneracy 3) and coexist with the gas phase, in 


f 

5 x 10“ 2 

6 x 10“ 4 (f op t) 

10 -5 

Ti 

0.50 

0.50 

0.50 

t 2 

1.35 

1.60 

1.82 

t 3 

1.70 

2.05 

2.33 

TABLE I. Temperature sets Tr= 3 , for BEG model, consider- 

ing the frequencies / obtained from free-energy weights [33]. 

f 

0.37 

0.11 (fopt) 

10" 4 

Ti 

0.50 

0.50 

0.50 

t 2 

1.12 

1.25 

1.70 

t 3 

1.24 

1.55 

2.14 

t 4 

1.34 

1.78 

3.20 

T 5 

1.42 

1.95 

4.34 

t 6 

1.50 

2.06 

6.10 


TABLE II. Temperature sets Tr= 6, for the BEG model, 
considering distinct frequencies / obtained from free-energy 
weights [331 ]. 


such a way that at p* the value po ~ 1/2 is verified. 
For this study, three values of R and frequencies were 
considered, whose temperature sets and distinct /’s are 
obtained from free-energy weights and are shown in Ta¬ 
bles HJ HU and [Hi] (BEG) and Tables EEI and ED (BL). 
Results are summarized in Figs. [Hand[2] for L = 20 and 
L = 18, respectively. As discussed previously, in all cases 
extremely low or large frequencies provide no precise re¬ 
sults, signaling trapping in a given phase (g = 0 or 1 for 
the BEG and p = 0 and 2/3 for the BL), due to low T’s 
or hardly exchanges. However, a similar result is verified 
for the intermediate / and the lowest R = 3, in such 
a way that only for intermediate frequencies (/ = f op t) 
when R = 6 and R = 8 , the crossing between coexist¬ 
ing phases is verified. However, for R = 6 the averages 
deviate greatly from its mean value 2/3, indicating that 
despite tunneling between coexisting phase occurs, it is 
not efficient. Improved results are obtained only for f opt 
and R = 8 , in which the tunneling is more frequent. 
These results present stark difference with those obtained 
by taking the free-energy weights (insets), in which for 
all replica numbers (including the lowest R = 3 case), 
it is possible to find an optimal / opt , ensuring proper 
tunneling between phases. Thus, the first evidence sug¬ 
gests that the choice of the replica number R is more 
important with approximated weights than with the free- 
energy ones. 

The results for the BL shown in Fig. [2] reinforce this 
idea. Only for the largest case R = 6 with / = f opt , some 
tunneling between the gas and LDL phases occurs. As 
for the BEG model, free-energy weights provide reliable 
results for all R' s (insets). The crucial importance of 
R is understood as follows: For two arbitrarily tempera¬ 
tures Ti and X) the exchange frequency obtained from the 
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f 

0.37 

0.21 (f opt ) 

10~ 4 

Ti 

0.50 

0.50 

0.50 

t 2 

1.12 

1.16 

1.70 

t 3 

1.24 

1.39 

2.14 

Ti 

1.34 

1.56 

3.20 

t 5 

1.42 

1.72 

4.34 

t 6 

1.50 

1.85 

6.10 

t 7 

1.58 

1.96 

6.60 

Ts 

1.64 

2.03 

7.00 


f 

0.37 

0.07 (f 0 pt) 

0.02 

Ti 

0.100 

0.10 

0.10 

t 2 

0.200 

0.23 

0.25 

t 3 

0.230 

0.30 

0.32 

Ti 

0.260 

0.32 

0.41 

Ts 

0.285 

0.39 

0.47 

t 6 

0.312 

0.43 

0.60 


TABLE VI. Temperature sets Tr=g, for the BL model, consid¬ 
ering distinct frequencies / obtained from free-energy weights 


TABLE III. Temperature sets Tr= 8, for the BEG model, 
considering distinct frequencies / obtained from free-energy 
weights 13311 . 


f 

0.025 

2 x 10” 4 (fopt) 

10 -5 

Ti 

0.10 

0.10 

0.10 

t 2 

0.25 

0.32 

0.35 

t 3 

0.33 

0.43 

0.48 


TABLE IV. Temperature sets Tr= 3 , for the BL model, consid¬ 
ering distinct frequencies / obtained from free-energy weights 


approximated weights is lower than those from the free- 
energy ones, whose difference becomes more pronounced 
for the optimal choices. For example, for the BEG and 
R = 3 the exchange between T\ = 0.50 and T 2 = 1.60 
(Table I) is performed with frequency f op t = 6.10~ 4 when 
one takes the free-energy weights. For the approximated 
weights, it reads 10 -6 . An efficient performance of the 
ST (e.g. the system visiting properly distinct coexist¬ 
ing phases) depends not only on frequent exchanges but 
also a reliable estimation of all temperature schedules, 
including the extreme temperature Tr. Thus, the com¬ 
promise between these points is reached (by taking the 
approximated weights) only for larger R' s. For exam¬ 
ple, in such regime (exemplified here for R = 8) the ex¬ 
change between Tf = 0.50 and T 2 = 1.16 (Table III) is 
performed with a considerable larger frequency, reading 
/ = 0.21 and / = 0.015 when one takes the free-energy 
and approximated weights, respectively. We remark that 
frequency 0.015 is larger than the value 10 -6 obtained be- 


f 

0.15 

0.01 (Ut) 

1.5 x 10" 3 

Ti 

0.10 

0.10 

0.10 

t 2 

0.20 

0.27 

0.29 

t 3 

0.25 

0.34 

0.39 

Ti 

0.29 

0.43 

0.50 


TABLE V. Temperature sets Tr= 4 , for the BL model, consid¬ 
ering distinct frequencies / obtained from free-energy weights 


R=3 R=6 R=8 




) -;-1-;- 0 - » , | "HI - U h I , - 

0 50 100 0 J§ 100 o 50 100 


t (lo 7 MCs) 


FIG. 1. Order parameter q as a function of t for the BEG 
model at T = 0.5 for L = 20. Distinct values of /’s and 
R’s were considered at the phase coexistence D* = 8.0000(1), 
with temperature sets obtained from Ref. 133) but using the 
approximated weights. The dashed lines correspond to steady 
value qo = 2/2>. The insets correspond to the results from 
free-energy weights. 


tween 0.5 and 1.60 (approximated weights). For the BEG 
model, in Fig. [3] we confirm these points by obtaining 
temperature schedules Tr's (for R = 6 and R = 8) with 
frequencies evaluated from the approximated wei ght s, in¬ 
stead of those available from free-energy values [33. As 
in previous cases, three frequencies have been considered, 
with temperature sets shown in Tables IVIII and IVIIII for 
R = 6 and R = 8, respectively. As in Fig. [I] larger /’s 
provide more frequent exchanges, but it is not efficient 
since the obtained Tr is low. On the contrary, lower / 
gives a larger Tr but temperatures exchanges are hardly 
performed (clearly shown for the lowest / = 0.021 when 
R = 6). By increasing R , it becomes possible to obtain 
an intermediate frequency f opt that fulfills the two re¬ 
quirements above for an efficient performance (larger Tr 
and frequent temperature exchanges). 

Fig. Ha) reinforces above ideas by showing for the 
BEG model the time decay of the order-parameter q 
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R=3 R=4 R=6 R=6 



R=8 

lF T 1 T - 

o t , , feQ - 22 - 



t(10 MCs) 



FIG. 2. For the BL model and T = 0.1, p vs t for the BL 
model at the phase coexistence p* = —1.6500(1) for distinct 
R’s and f’s, with temperature sets obtained from Ref. [3l| 
but using the approximated weights. Numerical simulations 
have been carried out for L = 18. The dashed lines denote 
the steady density po = 1/2 and the insets correspond to the 
free-energy weights. 


f 

0.18 

0.06 (/ opt ) 

0.02 

Ti 

0.50 

0.50 

0.50 

t 2 

1.02 

1.10 

1.15 

t 3 

1.31 

1.48 

1.58 

t 4 

1.53 

1.77 

1.88 

t 5 

1.68 

1.97 

2.08 

t 6 

1.83 

2.12 

2.47 


TABLE VII. Temperature sets Tr=6 , for the BEG model, 
considering distinct frequencies / obtained from the approxi¬ 
mated weights. 


FIG. 3. For the BEG model, T = 0.5 and L = 20, q 
versus t for distinct f’s and R’s at the phase coexistence 
D* = 8.0000(1), whose sets are obtained from approxi¬ 
mated weights. The dashed lines correspond to steady value 
qo = 2/3. 



f 

0.22 

0.10 (fo P t) 

0.06 

Ti 

0.50 

0.50 

0.50 

t 2 

1.00 

1.07 

1.10 

t 3 

1.27 

1.42 

1.48 

t 4 

1.47 

1.68 

1.77 

t 5 

1.63 

1.87 

1.97 

t 6 

1.77 

2.01 

2.12 

TV 

1.89 

2.15 

2.44 

T s 

1.98 

2.43 

3.03 


TABLE VIII. Temperature sets Th=s, for the BEG model, 
considering distinct frequencies / obtained from the approxi¬ 
mated weights. 


FIG. 4. For the BEG model, panel (a) shows the time decay of 
the order-parameter q for distinct f’s at the phase coexistence 
D* = 8.0000(1) for R = 8 and L = 20. In ( b ) q vs D for 
distinct L’s, with numerical results obtained from Tr=8 for 
fopt■ The continuous lines denote results for q obtained from 
Eq. ©. In (c) we compare results with those obtained from 
approximated (stars) and free-energy (circles) weights. 

starting from a fully occupied initial configuration for 
distinct f’s and R = 8. Note that the optimal choice 
fopt also ensures the faster convergence toward the steady 
value 2/3. In panel (6) we extend the obtained 7 h= 8 (for 
the coexistence point D*) for other D’s and distinct sys¬ 
tem sizes. As for free-energy case 0[i5> an unique Trs 
can be used for other control parameters and system sizes 
(panel (c)), providing to characterize very precisely the 
transition point, as predicted by Eq. (0) by taking only 
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FIG. 5. In (a) and (b), p versus t for distinct values of /’s and 
R’s at the phase coexistence point (T = 0.1, p* = —1.6500(1)) 
for the BL model and L = 18. The dashed lines denote the 
steady density po = 1/2. From top to bottom, intermediate 
temperature schedules read {0.183, 0.230}, {0.200, 0.266} and 
{0.30,0.40} (R = 3) and {0.183,0.230,0.273,0.312,0.348}, 
{0.200, 0.266, 0.322, 0.372, 0.412} and 

{0.300,0.400,0.445,0.526,0.670} (R = 6). In (c) we 
compare with results obtained from free-energy weights 
(circles). In ( d ) p versus p for distinct system sizes for 
T = 0.1. The continuous lines denote results for p obtained 
from Eq. (J3J). 


small system sizes. 

In Fig. [5] we compare distinct Tr’s for the BL model, 
by taking approximated weights. Confirming previous 
results, the optimal frequency f op t is obtained only by 
increasing R (panels (a) and (6)) with results equivalent 
with the free-energy weights (panel (c)). As for the BEG 
model, the same 7 r is extended to values of p and system 
sizes, as shown in panel ( d ) with results obeying Eq. ([3]). 

Next, we analyze the BC model. Besides testing the 
previous ideas, such analysis also aims to consider a 
phase transition ruled by the temperature. Thus, it 
also shows the reliability of the whole temperature set 
obtained from f opt . We take the value D* = 1.9968, 
in which results from the cluster algorithms Q, WL 
method @ and PT [l8| predicts a phase coexistence 
close to T = 0.4. Results for the tunneling between 
coexisting phases are considered for L — 16, as shown 
in Fig. [6] From top to bottom panels, tempera¬ 
ture schedules read {0.375,0.400,0.420,0.440,0.458}, 
{0.350,0.400,0.440,0.485,0.527} and 

{0.310,0.400,0.481,0.560,0.650} (R = 5) and 

{0.375,0.400,0.420,0.440,0.458,0.480,0.509,0.530}, 
{0.350,0.400,0.440,0.485,0.527,0.570,0.619,0.680} 
and {0.310,0.400,0.481,0.560,0.650,0.783,0.985,1.290} 
(R = 8). As in previous examples, small R leads to 
trapping of the system with the method performance 


FIG. 6. For the BC model, panels (a) and ( b ), q versus t for 
distinct values of /’s and R's at the phase coexistence point 
(T = 0.40, D* = 1.9968) for L = 16. The dashed lines denote 
the steady value qo = 2/3. In (c) we plot q vs T for distinct 
L’s. The curves cross at the transition point T = 0.40(1). In 
(d), we compare results obtained from the WL method and 
cluster algorithms. 


improving as R increases. As it can be seen, the 
achievement of an optimal 7 r ensures tunneling and 
faster convergence toward the equilibrium value. Also, 
the same 7 r is extended for distinct L’s and from Eq. 
© we see that all curves cross at T = 0.40(1) (panel 
(c)). For T > 0.5, we see an unusual behavior, signaling 
the large temperature regimes and the non-validity of 
Eq. ©. Comparison with cluster algorithms @] and 
WL results (32$ (panel (d)) confirms once again the main 
ideas that approximate weights can be used properly in 
order to give results equivalent with such tools. 

For the last example, Figs. [7] and [5] summarize results 
for the ALG model by taking the gas-LDL and LDL-HDL 
coexistence points for T = 0.30. In such cases, one ex¬ 
pects steady values close to po ~ 3/5 and <j >o « 0.425, re¬ 
spectively. The former value can be understood by recall¬ 
ing the LDL phase has density p = 3/4 and degeneracy 4 
and coexists with the gas phase. Since their weights are 
equal (1/5) at the phase coexistence, the steady po ~ 3/5 
is verified. The latter value is not easily understood, 
since the HDL phase is highly degenerated. Results for 
the tunneling between coexisting phases are shown for 
L = 12. As in all previous cases, extreme frequencies 
give results deviating from steady values, whereas the 
optimal frequency is signed for proper crossing between 
phases with accuracy improving as R increases. However, 
in such case each transition point study requires its own 
temperature schedule, since the gas-LDL coexistence line 
is shorter than the LDL-HDL, ending at respective dis¬ 
tinct critical temperatures T c = 0.55 and 0.825, respec¬ 
tively [29)]. For each 7 r we obtain the thermodynamic 
quantities, as shown in parts ©c) and ©(c). Note that 


































































FIG. 7. For the ALG model and L = 12, panels (a) and 
(b) show p versus t for distinct values of f’s and R’s at 
the gas-LDL phase coexistence point p* = —1.9986(1) and 
T = 0.30. The dashed lines denote the steady density po = 
3/5. From top to bottom, intermediate temperature sched¬ 
ules read {0.32, 0.34}, {0.406, 0.520} and {0.44, 0.56} (R = 3) 
and {0.32, 0.34, 0.355, 0.372}, {0.355, 0.395, 0.434, 0.473} and 
{0.38, 0.45, 0.53, 0.653} (R = 5). In (c) p versus p for distinct 
system sizes. The continuous lines denote results for p ob¬ 
tained from Eq. ©• In (d) we compare with results obtained 
from the PT (triangles) for L = 12. 


results are fully described by Eq. © (continuous lines). 
Also, both panels (d) show the equivalence of results 
with those obtained from the PT [35(. However, larger 
replica numbers and non adjacent replica exchanges were 
required for the system visiting properly the phases when 
the PT is considered [35j . 


V. DISCUSSION AND CONCLUSION 

Simulated tempering with approximated weights, pro¬ 
posed by Park et. al. [25|, presents a great advantage 
over other procedures, since the achievement of weights 
is readily obtained by executing simple and standard nu¬ 
merical simulations. Aimed at unveiling and optimiz¬ 
ing its performance, we scrutinized their main points 
focusing in the regime of first-order transitions at low 
temperatures, in which the existence of large trapping 
in metastable states makes the improvement of approx¬ 
imated weights a highly desired problem. Four systems 
having available results from parallel tempering, simu¬ 
lated tempering with free-energy weights, Wang-Landau 
method and cluster algorithms have been considered and 
thus they are important benchmarks for testing all ob¬ 
tained results. 

In all cases, results showed that approximated weights 


improves the efficiency of the ST as the replica number 
increases and optimized sets accelerating the convergence 
toward the steady regime and promoting frequent tun- 



FIG. 8. For the ALG model and L = 12, panels (a) and 
(6) show cj> versus t for distinct values of f’s and R’s at 
the LDL-HDL phase coexistence point p* = 2.000(1) and 
Ti = 0.30. The dashed lines denote the steady density 4>o = 
0.425. From top to bottom, intermediate temperature sched¬ 
ules read {0.55,0.70}, {0.60,0.78} and {0.65,0.84} (R = 3) 
and {0.50,0.579,0.650,0.705}, {0.55,0.70,0.795,0.840} and 
{0.60, 0.78, 0.89,1.20} (R = 5). In (c) <f> versus p for distinct 
system sizes. The continuous lines denote results for ob¬ 
tained from Eq. ©. In (d) we compare with results obtained 
from the PT (triangles) for L = 12. 


neling between coexisting phases have been found. We 
remark that the minimum value R ensuring proper sam¬ 
pling may depend on the system size and also the sys¬ 
tem temperature studied. We believe that the proposed 
recipe for the ST with approximated weights combined 
with Eq. © provides a very powerful method for deal¬ 
ing with discontinuous transitions, demanding relative 
short systems (as those considered here) and few control 
parameters [HI, [Hj. Although not necessary (in confor¬ 
mity with the main ideas concerning Eq. ©), the exten¬ 
sion for larger system sizes is straightforward. However, 
larger replica numbers or non-adjacent replica exchanges 
can be required. As a final comment, we remark that fu¬ 
ture method extensions include polymeric models and off- 
lattice systems, in which tempering methods have been 
extensively exploited, but the achievement of free-energy 
weights is difficult. These points should be addressed in 
future works. 
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